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Abstract 

Prospects are examined for the detection of a slow decay of the lightest neutralino 

(or any other longlived particles) at the CERN LHC and at Very Large Hadron 

Collider (VLHC). We first point out that such hadron colliders will become the 

“neutralino factory” producing 10®-10® neutralinos/yr, if gluinos and/or squarks 

actually exist below 0(1) TeV. The lightest neutralino (x?), usually assumed to be 

stable, will be unstable if lighter superparticles such as the gravitino (G) or axino 

(a) exist, or 7?-parity is not conserved. The decay signal would, however, be missed 

in usual collider experiments, particularly when the decay mostly occurs outside the 

detector. In order to search for such a slow decay of x5, we propose a dedicated 

experiment where the collision products are dumped by a thick shield, which is 

followed by a long decay tunnel. The decay product of x? can be detected by a 

detector located at the end of the tunnel. The slow arrival time and the large off 

angle (to the direction of the interaction point) of the decay product will provide 

a clear signature of slowly decaying Xi’s- One can explore the decay length (cr) in 

a wide range, he., 0.2 m to 1 x 10^ km for m^o = 25 GeV and 1 m to 2 km for 

m^o = 200 GeV at the LHC. This corresponds to the range of the SUSY breaking 

scale \Jf = 2 X 10® to 2 x 10^ GeV in case of the Xi ^ 7^5 decay predicted in 

gauge-mediated SUSY breaking models. At VLHC, one can extend the explorable 

range of m^o up to 1000 GeV, and that of pP up to ~ 1 x 10® GeV. In case of 
•^1 

the Xi 7® decay, the Peccei-Quinn symmetry breaking scale ijj can be explored 
up to ~ 5 X 10^^ GeV. The mass of the decaying particle can be determined by 
using the correlation between the energy and the arrival time of the decay product. 
With the setup we propose, one can also search for (i) other decay modes of x? such 
as R-parity violating one, (ii) slow decays of any other longlived neutral or charged 
particles, and (hi) heavy stable charged particles. 
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1 Introduction 


s = 14 TeV with luminosity C = 10^^ cm and Very Large 
s = 100-200 TeV with C = lO^^-lO^® cm-^s-^). 


The search for supersymmetric particles is now an integral part of all current, as well 
as future, experimental programs at high-energy colliders. Aside from many attractive 
features of supersymmetry (SUSY) [jl[, the driving force for these searches comes from 
the recognition that weak-scale SUSY, which is introduced to solve the gauge hierarchy 
problem, requires that the SUSY partners of the standard model (SM) particles must be 
accessible to experiments that probe the TeV energy scale. If this is the case, a large 
number of gluinos and squarks will be produced at future hadron colliders such as the 
CERN LHC (operated at 
Hadron Collider (VLHC; operated at 
Once produced, gluinos and squarks will subsequently decay to the lightest neutralino 
(x?). This means that such hadron colliders will become the “neutralino factory”, which 
is capable of producing up to 10® Xi’s per year. 

The Xi is usually assumed to be the lightest supersymmetric particle (LSP) and thus 
stable if i?-parity is conserved. It will then escape the detector, resulting in the famous 
missing energy signature for SUSY However, the Xi might not be altogether stable: 
If there exists another superparticle lighter than Xij such as the gravitino {G) or axino 
(a), the Xi will decay into, e.g., yG 0 or ya [§|. Such a light gravitino naturally exists 
in gauge-mediated SUSY breaking (GMSB) models [5-12] as well as in a class of no-scale 
supergravity (SUGRA) models [|^ . A light axino can also exist in SUGRA models with 
the Peccei-Quinn (PQ) symmetry [|^. As another possibility, the violation of R-parity 
leads to the unstable Xi even if it is the LSP. The Xi will then decay into such modes as 
qqv and yi/ [|^ . 


If the Xi decay takes place inside the detector, the resultant event topology would 
be very different from that in case of the stable Xv The experimental signatures have 
recently been studied for the y® decay into yG [16-25] and ya [^, motivated by the 
single eeyy -|- Ifj, event observed in the CDF experiment at the Tevatron . For the 
CDF event interpreted as the signal of the Xi decay, the inferred decay length cr is much 
below 1 m. However, subsequent reports on diphoton Ifj, distribution observed in the 
CDF and D0 experiments, as well as the analysis of the LEP data at \/s = 161 
GeV [0, do not give any further evidence for the y? decay into photon(s) with such short 
cr. Therefore, the possibility is still open for the Xi decay into yG or ya occurring mostly 
outside the detector. Actually, such a slow decay of Xi appears to be favored at least 
for relatively light Original GMSB models ]|^ prefer relatively high SUSY breaking 
scale, a/F > 10”^ GeV [0, |^, implying cr(x® —>• yG) > 100 km for m^o < 60 GeV. In 
case of the Xi 7® decay, the allowed range of the PQ symmetry breaking scale, 10® 
GeV < Ri < 10^^ GeV [^, |^, leads to cr > 10 km for m^o < 60 GeV. 

If the Xi decay is slow and mostly occurs outside the detector, the signature observed 
in usual collider experiments will be indistinguishable from that of the stable y®. Hence 
the signal of the Xi decay would be missed entirely. Even if a signihcant part of y^’s 
produced at hadron colliders decay into photon(s) inside the detector, it might be hard to 
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detect the signature with usual detectors, because a huge number of particles produced 
will make it difficult to identify a single photon not coming from the interaction point. 
In addition, the detectors for the approved experiments at the LHC are not designed to 
measure the direction of each photon precisely. 

Therefore, in order to search for a slow decay of Xi, we propose a dedicated experiment 
with the interaction point surrounded by a thick shield, by which most of the collision 
products are dumped. The y? would go through the shield and decay in a long decay 
tunnel. The decay product of y? can then be detected by a detector located at the end of 
the tunnel. We show that the slow arrival time and the large off angle (to the direction 
of the interaction point) of the decay product will provide an unambiguous signal of 
slowly decaying y^’s. We also show that, in case of the y? decay into a photon and a 
light invisible particle such as G and a, the mass of Xi can be determined by using the 
correlation between the energy and the arrival time of the decay photon. Furthermore, by 
reconstructing the decay kinematics with the estimated decay position, one can determine 
whether the hnal state is two-body or not. 

The remainder of this paper is organized as follows. We hrst show in Section ^ how the 
LHC can be the neutralino factory. In Section |], we briefly review the theoretical models 
of the Xi decay into a gravitino or an axino. Our strategy of the search for a slow decay 
of Xi is then described in Section |[ In Section |^, we illustrate the method of determining 
the mass of y? after the discovery of its decay. Section | is devoted to the discussion 
of the Xi decay search at VLHC. We also discuss in Section other searches which are 
possible with our setup, i.e., (i) other decay modes of Xi, (ii) slow decays of any other 
longlived neutral or charged particles, and (hi) heavy stable charged particles. Finally, 
we present our conclusions in Section ||. A realistic design for the detector is investigated 
with Monte Carlo simulations in Appendix In Appendix we estimate the rates of 
high-energy prompt neutrinos and muons, which can go through the iron shield and thus 
become the potential background. 


2 The LHC as the neutralino factory 


To simulate the production and decay at the LHC, we hrst generate various SUSY 


processes in pp collisions at ^/s = 14 TeV using ISA JET 7.13 with the CTEQ2L 
parton distribution functions 0]. In most cases, the dominant products are the gluino 
(g) and squark (g), which will then decay into y^’s through q ^ gq and g —»• qqXi if 
rriq > Trig, and g ^ qq and q qxi if If the gluino and squark are heavy ( > 1 

TeV), the production of gauginos (charginos and neutralinos) could become dominant. 
As an example, we take rriq = 1.45m^, p = —rrig and tan/3 = 2, and assume the gaugino 
mass unihcation, which is known to roughly hold not only in SUGRA models but also in 
most of GMSB models. 

Figure |I| shows the total cross sections for the production of gluinos, squarks and 
gauginos, as well as the anticipated numbers of events per 100 fb“^, at the LHC as a 
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function of rrig. We note that each SUSY event contains two Xi’s in the hnal state. Also 
shown in this hgure are the corresponding values of m^o ~ mg/T, which results from the 
gaugino mass unihcation. It is worth noting that the cross sections increase by a factor 
up to ~ 2 when the next-to-leading order QCD corrections are properly included . 


For the discussion below, we select three representative gluino and neutralino masses: 
nig = 300 GeV and m^o = 43 GeV (Gase 1); rrig = 550 GeV and m^o = 78 GeV (Gase 2); 
and nig = 800 GeV and m^o = 115 GeV (Gase 3). 

Figure || shows the angular distributions of x°’s produced at the LHG for (a) Gase 1, 
(b) Gase 2, and (c) Gase 3. The integrated luminosity is taken to be 100 fb“^, which 
corresponds to the runtime of 10^ s with C = 10^"^ cm“^s“^. The total numbers of 
produced 2 x 10®, 8 x 10® and 9 x 10® for Gase 1, 2 and 3, respectively. It can 

be seen that most of x®’s are produced in the forward direction. More events can thus be 
detected if one installs a detector at a smaller angle 6 with respect to the colliding proton 
beam. We here choose 6 = 25° (or cos6' = 0.906) for the central axis of the detector by 
compromising the rate of the signal and that of the potential background (high-energy 
prompt neutrinos and muons) as discussed in Appendix p. 

Figure shows the energy spectra of produced y^’s pointing to the detector which 
covers 6 = 25° ± 10° and an elevation angle of ±10°. We hnd that, in all three cases, 
the majority of y^’s are produced with relatively low energies, he., the energy spectra 
have a peak at E^o/m^o ~ 1.5 to 2. This means that most of the produced x5’s are not 
extremely relativistic, having (3 = v/c < 0.96. This makes the time measurement of the 
decay products a powerful tool for the identihcation of slowly decaying ;y®’s. 


3 Models of the lightest neutralino decay 

3.1 The decay into a gravitino 

Supersymmetry (SUSY) is widely acknowledged to be the best motivated extension of 
the SM. This is not only because it provides a natural solution to the gauge hierarchy 
problem, but also there exist many pieces of indirect evidence, such as the correct pre¬ 
diction of sin^6*iy, the convergence of the SM gauge couplings at very high energies, and 
the heaviness of the top quark as required by the radiative breaking of the electroweak 
gauge symmetry. Moreover, SUSY can provide a plausible candidate for cosmological 
dark matter. These hints, as a whole, appear to constitute circumstantial evidence for 
SUSY. 

SUSY, however, must be broken because no SUSY partners of ordinary particles have 
been discovered. The breaking of SUSY is usually conceived to occur in a hidden sector 
at some very high energies, which will probably remain unaccessible to direct collider 
experiments at least in the near future. Thus the investigation into the origin of the 
SUSY breaking is of great importance since it could be the unique window to look into 
the world of such high energies, perhaps up to the Planck scale. Two scenarios have 
been proposed on how the SUSY breaking is communicated from the hidden sector to 
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the observable sector: In one scenario, as in supergravity (SUGRA) models |Q, SUSY 


breaking is mediated by gravity. In this case, the SUSY breaking scale \/F is necessarily 
of order 10^^ GeV to give the superparticle masses around the weak scale. The other is 
the gauge-mediated SUSY breaking (GMSB) scenario [5-12], where SUSY is broken at a 
scale as low as 10^-10^ GeV, and the gauge interactions act as the messengers of SUSY 
breaking. Because of this relatively low value of \/F, the mass degeneracy among squarks 
or sleptons, which results from the fact that gauge interactions are flavor blind, is hardly 
broken by the evolution to the weak scale, in contrast to SUGRA models. This ensures 
sufficient suppression of ffavor changing neutral current (FGNG) processes. In addition, 
GMSB models may also provide a solution to the SUSY CP problem p. 

These two scenarios also lead to distinctive consequences concerning cosmological dark 
matter. In SUGRA models, the lightest neutralino (y?) is usually the LSP, and can be 
the dominant component of dark matter. On the other hand, GMSB models naturally 
predict that the gravitino {G) comes out to be the LSP. The gravitino with mass ~ 0.5 

iHj. We note that it 


keV is claimed to be a good candidate for warm dark matter 
is practically impossible to directly detect the gravitino dark matter in laboratories, in 
contrast to the neutralino dark matter, whose direct detection is being pursued in many 
experiments. 

In this way, for hnding the future direction of particle physics and cosmology, it is of 
crucial importance to determine which scenario is realized in nature. As will be shown 
later, our experiment to search for a slow decay of Xi can probe a wide range of cr, which 
corresponds to y/P of order 10® to 10® GeV in case of the Xi iG decay. This covers 
much of the y/P range presently considered in GMSB models. 

The gravitino mass ttiq is related to y/P (on the condition of the vanishing cosmological 
constant) by 

P ^ ( y/F V 


= 


■y/SAfp 


= 2.13 


3 X 10® GeV 


keV, 


( 1 ) 


where Mp = = 2.44 x 10^® GeV is the reduced Planck mass. This leads to 

niQ ~ 0.24 and 24 keV for y/P = 10® and 10^ GeV, respectively. Note that the gravitino 
heavier than 0(1) keV overdoses the Universe in the standard cosmology. Such a heavy 
gravitino, however, is still allowed to be the dark matter if the reheating temperature 
after the inflation is low enough ||, . 


The longitudinal component of the gravitino, he., the Goldstino, couples to matter 
with weak (not gravitational) interaction strength proportional to P~^. Through this 
coupling, the Xi will decay to a photon and a gravitino with a partial width of 0 


r(x? - xG) = 


_A-l 

167rF2 ’ 


( 2 ) 


where is the square of the photino admixture of the Xi- If fhc Xi is a pure bino, 
Ki^ = cos^ 6w — 0.77. This partial width is most likely to be the dominant one if the Xi 
is the next-to-lightest supersymmetric particle (NLSP) with not too small ki^, resulting 
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in the decay length 
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3 X 106 GeV 
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( 3 ) 


For > 3 X 10® GeV and m^o <100 GeV, this decay length is much longer than a 
typical detector size, leading that most of the decays will occur outside the detector. It 
has recently been pointed out |jl|, that such large a/F ( > 10^ GeV) is expected in 
original GMSB models p. In addition, other GMSB models [9-11] tend to predict that 
SUSY breaking occurs at relatively high ( > 10® GeV) energies, although y/F = 10®-10® 
GeV can naturally £t in with a very recent model |^. Equation @ also indicates that 
the measurement of both m^o and cr will give direct information on y/F. 


3.2 The decay into an axino 

The Peccei-Quinn (PQ) symmetry was introduced to solve the strong CP problem |40| . 
This symmetry should be explicitly broken, leading to a light pseudoscalar, the axion (a), 
whose mass is related to the PQ symmetry breaking scale F^ by 


m 


a 


6.2 X 10“^ 


Fa/N 

109 GeV 


-1 


eV, 


( 4 ) 


where N is the QGD anomaly factor of the PQ symmetry. The axion is one of good 
candidates for dark matter owing to its feeble interaction with matter. The allowed range 
of Fa is conhned to 109-10^^ GeV Q from laboratory experiments as well as from 
astrophysical and cosmological considerations. 

If both the PQ symmetry and SUSY hold in nature, there should exist the fermionic 
superpartner of the axion, the axino (a). The axino can remain light even after the SUSY 
breaking in SUGRA models fl^. In particular, the axino with mass ~ 0(1) keV could 


serve as warm dark matter as is the case for the gravitino. 

If such a light axino exists, the Xi can decay into ya. In case of the Xi being a pure 
bino, the decay width can be written as @ 


r(x? ^ 7a) = 

which corresponds to the decay length 

CT = 3.6 


25a^ 


m?o 

Xi 


11527r3cos2 Ow {Fa/N) 


2 ’ 


( 5 ) 


Al 


100 GeV 


-3 


' Fa/N ' 
109 GeV 


km. 


( 6 ) 


The lower limit on the PQ symmetry breaking scale. Fa > 10® GeV, implies cr > 10 km 
for m^o < 60 GeV. 

Xi ~ 


5 


























4 Strategy of the neutralino decay search 

4.1 Dedicated experiment 

We propose a completely new type of collider experiment, which is schematically shown in 
Fig. Gluinos and squarks (and possibly gauginos) are produced in pp collisions at the 
interaction point and promptly decay into x5’s. The y^’s then enter the decay tunnel after 
going through the shield, by which most of the collision products are dumped. As a simple 
example, we here consider an iron shield of 10 m thickness, and take the tunnel length to 
be L = 41.6 m, although a much longer tunnel is certainly possible. The acceptance of 
the tunnel is taken to be 0 = 25° ± 10° and an elevation angle of ±10°. We choose this 
angle 9 by compromising the rate of the signal and that of the potential background (see 
Appendix P). Two layers of the anti-coincidence hodoscopes possibly with track detectors 
are installed just on the back end of the shield for the rejection of punchthrough muons. 
It can also be used for the trigger of heavy longlived charged particles. In addition, the 
scintillation counters and/or track detectors can be installed on the wall around the decay 
tunnel to reject potential background events produced by cosmic rays. If the Xi decays to 
a photon and a light invisible particle in the tunnel, the decay photon will reach the end 
of the tunnel and enter a shower detector which has a spherical front face, and initiate 
the electromagnetic shower. The incident position, angle, and energy of the photon, as 
well as the arrival time relative to the RF clock of the accelerator, are measured by the 
detector. 

Instead of the bulk iron shield, we can use (at extra costs) an “active” dvr-shield, 
consisting of segmented calorimeters, magnetized iron plates and the muon trackers. This 
setup will enable us to measure the total (and missing) transverse energy of the event and 
to tag the existence of high-energy muons. 

4.2 Event simulation and the discovery signal 

The decay probability of Xi within a distance of x is given by 

P{x) = 1 - (7) 

where 7 = E^o/m^o is the Lorentz factor, and cr is the decay length. As a typical case, 
the probability that the with 7 = 2 decays in the tunnel (10 m < a: < 51.6 m) is 
calculated to be 20% for cr = 100 m, and 2.4% for cr = 1000 m. 

In the following Monte Carlo simulations, we assume that the decay photon is emitted 
isotropically with energy equal to m^o/2 in the rest frame of the Xi^ is the case for the 
Xi —^ 7 G decay with itIq -C photon 4-momentum p^) in the laboratory 

frame is then calculated by the Lorentz boost. We show in Fig. ^ a typical generated 
event of the decay into a photon and a light invisible particle. 

Figure shows the distributions of the Lorentz factor 7 for y^’s entering the decay 
tunnel for Case 1, 2 and 3 in (a), (b) and (c), respectively. The numbers of Xi’s entering 
the tunnel are 9 x 10®, 4 x 10® and 5 x 10^ for Case 1, 2 and 3, respectively, for total 
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integrated luminosity of 100 fb“^, which we assume for all the simulation plots shown in 
this paper. The distributions of 7 for decaying in the tunnel are also shown in Fig. |^. 
We here take as an example the decay lengths calculated by Eq. (|^) for a/F = 3 x 10® 
GeV, he., cr = 7.0 km, 350 m, and 52 m for m^o = 43 GeV (Gase 1), 78 GeV (Gase 2) 
and 115 GeV (Gase 3), respectively. The resultant numbers of x®’s decaying in the tunnel 
are 1.4 x 10^, 1.3 x 10^ and 8 x 10^ for Gase 1, 2 and 3, respectively. Also shown in this 
hgure by hatched histograms are the distributions for x®’s whose decay photons enter the 
shower detector. The corresponding energy spectra of the decay photons are shown in 
Fig.0 for (a) Gase 1, (b) Gase 2, and (c) Gase 3. The spectra have a peak at ~ m^o 
and extend up to ~ 1000 GeV. 

The clear separation of the Xi decay signal from the potential background, he., high- 
energy prompt neutrinos and muons as well as the cosmic-ray events, can be attained by 
measuring the arrival time and the off angle ip to the direction of the interaction point (see 
Fig.|). Figure | shows the scatter plots of ip vs. the arrival time for (a) Gase 1, (b) Gase 2, 
(c) Gase 3, and (d) background neutrinos and muons (i^/p). The background z//p events 
all have ip ~ 0°, and arrive in narrow bunch structure (here we assume at = 0.2 ns). On 
the other hand, most of the decay photons have large ip (up to ~ 30°) and arrive with 
signihcant delay (typically 1-20 ns), because most of the produced Xi’s are not extremely 
relativistic. We show in Fig. ^ the arrival time distributions. The evidence for a slow 
decay of the heavy parent particles is clear in the peak position, which is displaced by 
about 2 ns from the exact bunch position, as well as in the long tail after the peak. 

The decay length cr can in principle be estimated from the distribution of the events 
in the scatter plot of ip vs. the arrival time. Moreover, one could measure cr from the 
variation of the signal count rate correlated with the change of L, which can be performed 
by moving the detector on rails. Together with the measured value of m^o (see Section 
one can determine the strength of the interaction which governs the decay. In partic¬ 
ular, if one interprets the observed events as the —>• 7G (y? —>• 7a) decay, the SUSY 

breaking scale \/F (the PQ symmetry breaking scale Fa) can be determined. 

4.3 Explorable range of the decay length at the LHC 

The large off angle ip and the slow arrival time are the generic features of slowly decaying 
heavy particles, and will thus always provide a clear signature of slowly decaying x®’s. 
Based on this prospect, we now estimate the explorable range of the decay length cr. 

Figures and ^ represent typical examples of the off angle ip and the arrival time 
distributions for cr much longer than (Gases 1 and 2) or comparable to (Gase 3) the 
tunnel length. For longer cr, the characteristic shape of the distribution will not change, 
while the number of detected events will decrease in proportion to (cr)“^. In order 
to estimate the longest explorable decay length, we assume that 10 events with such a 
characteristic distribution are enough to discover the signal. 

For cr much shorter than the tunnel length, most of the detected events will be due 
to the Xi decaying immediately after going through the shield. Such events will show the 
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smaller off angle and shorter arrival time as compared to the cases shown in Figs. ^ and 
p. The off angle will nevertheless have hnite non-zero values (typically of 1 to 10°), and 
the arrival time will dehnitely be deviated (by 1 to 10 ns) from that of background v/^ 
events, because the y? flies at least the shield thickness before decaying. In estimating 
the shortest explorable decay length, we assume 100 such events are enough to discover 
the signal. 

Figure shows the resultant explorable range of cr for various m^o at the LHC with 
300 fb“^. It can be seen that the explorable range of cr is wide, i.e., 0.2 m to 1 x 10 ® 
km for m^o = 25 GeV, and 1 m to 2 km for m^o = 200 GeV. We also show in this hgure 
the predicted curves of cr for the Xi decay into yG and jd for typical values of y/F and 
Fa, respectively. We hnd that one can explore y/F in a wide range of 2 x 10 ® to 2 x 10 ^ 
GeV through the Xi 7^ decay. This range includes an interesting case of y/F ~ 10 ® 
GeV, where the gravitino has mass ~ 0.5 keV and can be the dominant component of 
dark matter in the standard cosmology. In case of the Xi 7® decay. Fa can be explored 
up to ~ 10 ^® GeV for m^o < 60 GeV. 


5 Determination of the neutralino mass and decay 
kinematics 


After discovering the signal of the Xi decay, one can determine the mass of Xi by using 
the correlation between the energy and the arrival time of the decay photon, if enough 
number of decay events can be accumulated. 

For this purpose, we hrst extract the approximate time-of-flight (TOF) of Xi, ^tof, 
from the time difference between the signal of the shower detector and the RF clock of the 
accelerator, as illustrated in Fig. n. In this process, we assume that each Xi comes from 
the closest possible bunch of the proton beam on the condition of /? < 1. The probability 
of attributing an event to a wrong starting bunch is small, i.e., 20-25% for the three cases 
described above. Using the extracted value of tTOF, we dehne the parameter (3^ as 


[ 3 * = 


I OF I 

C^tOF’ 


( 8 ) 


where | OF | is the distance from the interaction point to the front face of the shower 
detector. The distributions of (3~^ are shown in Fig. ^ for (a) Gase 1, (b) Gase 2, 
(c) Gase 3, and (d) background vj There is an upper cutoff at I3~^ ~ 1.15, which 
results from the modulation of ^tof by the bunch spacing of 25 ns. The background v/^ 
events make a sharp peak at j3~^ ~ 1 with no tail in contrast to the decay signal. 

Figure p!3| shows the scatter plots of E.y vs. 7* = ( 1 -/ 3 ^)“^/^ for (a) Gase 1 , (b) Gase 2, 
and (c) Gase 3. The photon energy resolution is assumed to be cx^/F = 33%/ ^Je (GeV). 
The lower cutoff of 7* at ~ 2 corresponds to the above-mentioned cutoff of (3~^ at ~ 1 . 15 . 
One can also hnd an edge with constant ratio of E.y to 7*. This edge is composed of the 
events in which the decay occurs near the front face of the shower detector with the decay 





photon emitted toward the flight direction of the Xv For such events, | OE | and tTOF are 
nearly the same as the actual flight length and TOF of the Xi, respectively, and the decay 
photon has the full energy of the Xi^ -^7 ~ resulting in the edge corresponding 
with From the position of this edge, the mass of Xi can be determined 

with an accuracy of < 20%. For this purpose, the shower detector should be designed 
to measure the photon energy in a wide range of 10 to 1000 GeV. Figure |T^ shows the 
corresponding distributions of E^/x* for Cases 1 to 3. An upper endpoint at E^/x* = m^o 
can be seen in each case. 

We now try to determine the nature of the Xi decay, he., whether it is two-body or 
not. Although the symmetric peak at E^/x* = m^j2 in Fig. |T^ may be the hrst hint 
of the two-body decay, we can prove the two-body nature by directly reconstructing the 
decay kinematics. For this purpose, we hrst select events with i\) > 10° based on the ifj 
distribution shown in Fig. |T^ for (a) Case 1, (b) Case 2, (c) Case 3, and (d) background 
z///i. After this V’ cut, only events in which the Xi decays relatively near the shower 
detector can survive. The efficiency of this cut is ~ 40% for each of the three cases. 

Figure shows the scatter plots of E^ vs. 7* for Cases 1 to 3 after the cut oi'ip > 10°. 
It can be seen that only photons with relatively low energies survive. Thus the two 
methods of determining the Xi mass, he., the use of the edge in the scatter plot of E^ 
vs. 7*, and the direct reconstruction of the decay kinematics, are complementary in that 
different events are mainly used in these two methods. 

We then simply assume that the decay position D is the middle point between C and 
E, both of which are dehned in Fig. Using this assumption, the velocity of the Xi can 
be approximately expressed by 

/?obs =-, (9) 

c^xoF — I DE I 

where | CD | and | DE | represent the flight path of the Xi and that of the photon, 
respectively. One can then calculate the reconstructed mass as 

^7G = ^7(1 - f^ohs COS0iab)(l - (10) 

where ^lab is the angle between the momentum of the Xi and that of the photon (see Fig. 
Figure 0 shows the distributions of for (a) Case 1, (b) Case 2, and (c) Case 3. It 
can be seen that the Xi mass can be determined with an accuracy of ~ 30% for each case. 
If one hnds a peak in the distribution of at the right position expected from the edge 
in the scatter plot of E^ vs. 7*, one can conclude that the observed events are really the 
signature of the two-body decay of Xi mto a photon and an invisible particle. 

6 Prospects for the neutralino decay search at VLHC 

Very Large Hadron Collider (VLHC) has been discussed recently as a post-LHC 
machine, being operated at ^/s = 100-200 TeV with C = lO^^-lO^® cm“^s“^. We stress 
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here that the VLHC will become the neutralino factory much more powerful than the 
LHC. Figure ^ shows the total cross sections for the production of gluinos, squarks and 
gauginos, as well as the anticipated numbers of events per 1000 fb“^, in pp collisions at 
a/s = 100 TeV as a function of rrig. Also shown in this hgure are the corresponding values 
of m^o ^ mg/T. 

Figure shows the explorable range of cr for the Xi decay at the VLHC with 3000 
fb“h As in Fig. 13. the upper end of the explorable region corresponds to the level of 
10 decay events detected, and the lower one is derived from the condition that 100 decay 
events are detected. The explorable range of cr is found to be 0.1 m to 1 x 10^ km for 
m^o = 25 GeV, 0.2 m to 5 x 10^ km for m^o = 200 GeV, and 1 m to 1 km for m^o = 1000 
GeV. We also show in this hgure the predicted curves of cr for the Xi jG and Xi 7® 
decays for typical values of VF and respectively. It is found that one can extend the 
reach for y/F up to 1 x 10® GeV through the Xi iG decay. In case of the Xi 7® decay, 
Fa can be explored up to ~ 5 x 10^^ GeV for relatively light Xv If is worth noting that 
our experiment is unique in that one can explore such a high value of Fa in laboratories. 
In addition, our experiment can be performed under very high luminosity exceeding 10®® 
cm“^s“®, in contrast to usual collider experiments. 


7 Other possible searches 


7.1 Other decay modes of the lightest neutralino 

One can search for other modes of the Xi decay by using the setup proposed in this paper. 
If the x\ is sufficiently heavy, it can decay into Z^G or h^G (where is the Higgs boson) 
in the gravitino LSP scenario [^, In case the Xi is higgsino-like, the h^G decay can 
even be dominant. Since the and bosons subsequently decay into qq or , the 
visible signature of these Xi decays will be two hadron jets or two leptons without any 
signal in the anti-coincidence hodoscopes. With the help of track detectors at the end of 
the decay tunnel, the decay point can be exactly located, which makes the kinematical 
analysis of the decay and the mass determination much more streightforward than in case 
of the Xi iG decay described in Section |^. The Xi can also decay into such modes 
as qq'l^, qqv and yz/ through weak violation of F-parity [O. Furthermore, the decay 
of other longlived neutral particles, e.g. heavy neutrinos, can be detected with the same 
setup and method. 


7.2 Slow decays of longlived charged particles 

Our setup can also be used to search for slow decays of heavy longlived charged parti¬ 
cles. Although the existence of such particles would in principle be indicated in collider 
experiments by charged particle tracks with ionization larger than the minimum one, the 
signature of their decay will be missed if it occurs mostly outside the detector. On the 
other hand, with our setup, one can clearly detect the decay signal of such longlived par- 
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tides as in case of the Xi decay. The “anti-coincidence” hodoscopes can be used as the 
trigger counter. 

There are several candidates for longlived charged particles in SUSY models. For 
example, if the gravitino or axino is the LSP, a charged superparticle such as the slepton 
(/“) can be the NLSP, and may decay into l~G with cr much longer than the detector size. 
The chargino (xr) can be another example of the longlived NLSP if the mass difference 
between the xL cind the x? (assumed to be the LSP) is so small that it can decay only via 
xr e~UeXi [H]- Weak violation of i?-parity would also allow the existence of longlived 
charged particles |Q. In this case, the slepton can be the LSP, decaying slowly via, e.g., 
l~ l~v. Furthermore, heavy longlived charged particles could also exist in the SUSY 
breaking sector and messenger sectors of GMSB models o, There are also many 
other speculations on heavy longlived charged particles [46-49]. 


7.3 Heavy stable charged particles 

If the cr of heavy charged particles is extremely long, they will appear to be “stable” 
particles even in our experiment. One can search for such particles with our setup by 
precisely measuring their time-of-flight over 40 m distance between the “anti-coincidence” 
hodoscopes and the shower detector at the end of the tunnel. The j3~^ distribution for 
heavy stable charged particles is expected to have a long tail following the peak at ~ 1, 
in contrast to muons which make a sharp peak at = 1. Long flight path in our setup 
leads to excellent resolution of /3, which will provide a clear signal for the heavy stable 
charged particle. One could then determine its mass by measuring its momentum with 
magnetized iron and tracking chambers, which can be installed either in the active shield 
or at the very end of the shower detector. 


8 Conclusions 

We have investigated the detection of a slow decay of the lightest neutralino (or any other 
longlived particles) at the CERN LHC and at VLHC. After pointing out that such hadron 
colliders will become the “neutralino factory” producing 10®-10® Xi’s/yr if SUSY particles 
actually exist below 0(1) TeV, we have shown that a slow decay of x? can be detected in 
a dedicated experiment in which the collision products are dumped by a thick shield. The 
decay product of x? can then be detected by a detector located at the end of a long decay 
tunnel. The slow arrival time and the large off angle V’ of the decay product will provide 
a clear signature of slowly decaying Xi’s. Considering the Xi decay into a photon and a 
light invisible particle such as the gravitino (G) or axino (a), we hnd that one can explore 
the decay length (cr) of x? in a wide range, he., 0.2 m to 1 x 10^ km for m^o = 25 CeV 
and 1 m to 2 km for m^o = 200 CeV for the case of the LHC. This corresponds to the 
range of the SUSY breaking scale y/F = 2 x 10® to 2 x lO’^ CeV in case of the x? —^ xG 
decay. At VLHC, one can extend the explorable range of m^o up to ~ 1000 CeV, and 
that of \/F up to ~ 1 X 10® CeV. In case of the Xi 7® decay, one can explore the PQ 
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symmetry breaking scale up to ~ 5 x 10^^ GeV. The mass of the decaying particle 
can be determined with an accuracy of < 20% using the correlation between the energy 
and the arrival time of the decay photon. With the setup we propose, one can also search 
for (i) other decay modes of Xi such as i?-parity violating one, (ii) the decay of any other 
longlived neutral or charged particles, and (hi) heavy stable charged particles. 

The experiment we proposed in this paper is a completely new type of collider ex¬ 
periment, which may turn out to be of crucial importance. The civil engineering of the 
intersection region at future hadron colliders has to be designed so that it can accommo¬ 
date long decay tunnels for such experiments. 
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Appendixes 

A Realistic design for the shower detector 

Figure ^ shows a schematical view of a realistic design for the shower detector which 
we have investigated. To separate electromagnetic showers initiated by photons from the 
background, we design the detector consisting of the preshower detector (PS), two parts 
of the electromagnetic calorimeter (EMI and EM2), three planes of plastic scintillators 
(SI, S2 and S3), the lead absorber (LI), and the hadron calorimeter (HAD). In addition, 
we can install the tracking chambers both in front of and behind the whole detector for 
the dehnite detection of high-energy muons (and possibly heavy stable charged particles). 

The preshower detector (PS) comprises two planes of multiwire proportional cham¬ 
bers (MWPCs) with cathode pad readout, with 2Xq and IXq lead converter in front, 
respectively (radiation length Xq being 0.56 cm for lead). The scintillator SI in front of 
PS can work as the anti-coincidence counter for charged particles. The scintillator S2, 
located just behind PS, is used for triggering. The electromagnetic calorimeter is divided 
into two parts, each of which comprises modules with 10 x 10 cm^ transverse dimen¬ 
sions. The hrst part (EMI), which is located at 170 cm downstream of PS, consists of 5 
layers of lead (IXq) and scintillators (1 cm). The second part (EM2), which is located 
at 220 cm downstream of EMI, consists of 20 layers of lead (IXq) and scintillators (1 
cm). Total radiation lengths of EMI and EM2 are 5.1Xo and 20.4Xo, respectively. The 
thick scintillator S3, located between EMI and EM2, is used not only for triggering but 
also for the measurement of the arrival time. The hadron calorimeter (HAD) contains 
15 layers of 5 cm thick iron and 2 cm thick scintillators. Total thickness is ~ 5 nuclear 
interaction lengths. Hadronic showers will deposit much of its energy in HAD. Between 
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EM2 and HAD, the lead absorber LI of 20Xo in thickness is inserted to totally absorb 
electromagnetic showers. 

A coincidence of two planes of scintillators, i.e. S2 and S3, constitntes the hrst trigger. 
We also reqnire the signals from EMI and EM2 to exceed a threshold, which leads to 
the rejection of (i) mnons going throngh the detector as minimnm ionizing particles, and 
(ii) hadronic showers initiated by high-energy nentrinos in HAD. 

For this detector conhgnration, we have performed Monte Carlo simnlations nsing 
GEANT 3.21 a. A typical event for 100 GeV photon incident on the shower detector 
is shown in Fig. |^. We can see that the electromagnetic shower starts to evolve in PS, 
reaches its maximnm after passing EMI, and is mostly absorbed in EM2. Fignre p2| shows 
the distribntions of energy deposited in active layers of (a) S2, (b) EMI, (c) EM2, and 
(d) HAD, by 100 GeV photons and charged pions which are incident at 0° on the shower 
detector. It can be seen that electromagnetic showers (initiated by photons) can clearly 
be separated from hadronic showers (initiated by pions) nsing the deposit energy of each 
component. 

The energy of incident photons can be determined from the snm of deposit energy in 
active layers of EMI and EM2, denoted as Eemi + .Eem 2, whose distribntions are shown 
in Fig. ^(a) for photons incident at 0° with energy E = 20, 50, 100 and 200 GeV. The 
inefficiency concerning the scintillation light collection is not inclnded. Fignre |^(b) shows 
the relation between E and E^mi + -Sem 2 for 0° and 30° incidence. It can be seen that 
good linearity is realized for E = 20-200 GeV. The resnltant energy resolntion of the 
shower detector as a fnnction of incident photon energy is shown in Fig. ^(c). The fit 
resnlt for 0° data is (Je/E = 32.6 %/E (GeV) for E = 20-200 GeV. 

The shower direction is obtained from the reconstrncted position of the shower in PS 
and that in EM2. The position of the shower in PS can be determined from the pad 
signal of MWPGs, whose spatial resolntion is assnmed to be 5 mm. The shower center 
in EM2 can be estimated from the observed division of energy within a local cluster of 
modules. We hnd that good spatial resolution ~ 2 cm) in comparison with the size 
of the module (10 x 10 cm^) can be obtained owing to the large gap between EMI and 
EM2, which makes the shower spread widely before entering EM2. The resultant angular 
resolution averaged over various incident positions is shown in Fig. |2^. It is found to be 
ae < 0.3° for photons with E = 50-200 GeV. The angular resolution could be improved 
if tracking chambers are used as active layers of EMI and EM2 instead of scintillators. 


B Rates of high-energy neutrinos and muons 

Although a large number of particles will be produced at the LHG, most of them (hadrons, 
and 7) will be dumped in our setup by the iron shield with thickness of ~ 60 nuclear 
interaction lengths. However, high-energy neutrinos and muons can go through the shield 
and might interact with the material of the detector. The main source of high-energy 
neutrinos and muons is expected to be “prompt”, i.e., come from semileptonic decays 
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of bottom and charm hadrons |^, while the contribntion from the decay of vr and K 
mesons is negligible at least in onr setnp. In order to estimate the rates of these prompt 
nentrinos and mnons, we have performed Monte Carlo simnlations using ISA JET 7.13 
3B| for pp collisions at y/s = 14 TeV, by taking the total cross section of the bottom plus 


charm production to be 2 mb. Total 2.4 x 10® events are generated, and the neutrinos and 
muons which point to the detector covering a polar angle of ^axis ± 10° and an elevation 
angle of ±10° are recorded. We consider four cases, he., 6'axis = 15°, 20°, 25° and 30°. 


B.l Interacting neutrinos 

The energy spectra of the prompt neutrinos are shown in Fig. The spectra have a 
long high-energy tail when ^axis is small. The neutrino could mimic the yG or ya signal 
if it interacts in the preshower detector and deposits enough energy (see Appendix ^ for 
the detailed design for the shower detector). As a simple estimate of such background, 
we calculate the number of neutrinos which have energy > 20 GeV and interact with 
3Xo ~ 19 g/cm^ lead. The interaction cross section is taken to be 5 x 10“^®(Ej,/GeV) 
cm^/nucleon. The estimated numbers are ~ 0.2, ~ 0.02, < 0.01 and < 0.001 per 100 
fb“^ for 6'axis = 15°, 20°, 25° and 30°, respectively. Thus the background induced by the 
prompt neutrinos is negligibly small. 


B.2 Muon-induced background 


Figure ^ shows the energy spectra of the prompt muons. Only the muons with > 20 
GeV can go through the shield because of the ionization energy loss. The resultant rates 
of muons entering the tunnel are then estimated to be ~ 5 kHz, ~ 600 Hz, < 200 Hz and 
< 80 Hz for 6*axis = 15°, 20°, 25° and 30°, respectively, with luminosity of 10^^ cm“^s“^. 
The single rate of the prompt muons will thus be comfortably small if we take ^axis ^ 25°. 


High-energy muons can deposit a signihcant portion of their energy by bremsstrahlung 
or direct e’''e“ pair production, both of which could lead to the electromagnetic shower. 
When combined with inefficiencies in the anti-coincidence counters, these muons might 
in principle mimic the yG or ya signal. Folding the simulated spectra by the showering 
probability of muons, we estimate the numbers of muon-induced high-energy {E > 10 
GeV) showers starting to evolve in the preshower detector to be ~ 10^, ~ 10®, < 2 x 10® 
and < 10® per 100 fb“^ for 6*axis = 15°, 20°, 25° and 30°, respectively. Assuming a 
veto inefficiency of ~ 10“® by the staggered layers of anti-coincidence scintillators, the 
muon-induced background is expected to be negligibly small for ^axis ^ 25°. 


Even if the rate of such muon-induced background turns out to be non-negligible, 
off-line analysis will eliminate such background based on the information on the arrival 
time as well as on the off angle 'ijj. 
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Figure 1: Total cross sections for the production of gluinos, squarks and gauginos in 
pp collisions at ^/s = 14 TeV as a function of the gluino mass rrig. The corresponding 
mass of the lightest neutralino m^o is also shown. Also given in the vertical scale is the 
corresponding number of events per 100 fb“^. 
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Figure 2: Angular distributions of Xi’s produced at the LHC with 100 fb ^ for (a) Case 1, 
(b) Case 2, and (c) Case 3. 


19 










E .0 (GeV) 

Xi 


Figure 3: The energy spectra of which point to the detector for (a) Case 1, (b) Case 2, 
and (c) Case 3. 
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Figure 4: Schematic view of the setup. The lightest neutralino (x?) is produced in pp 
collisions, goes through the iron shield, and decays in the tunnel to a photon (7) plus 
a light invisible particle. The decay photon will then enter the shower detector, and be 
detected. An “active” 47r-shield consisting of segmented calorimeters, magnetized iron 
plates and the muon trackers can also be used instead of the bulk iron shield. 


21 





































































































































































































































Figure 5: An example of the simulated Xi decay into a photon and a light invisible 
particle. Dash-dotted, dotted and solid lines indicate tracks of the y?, > 50 MeV photons 
and electrons, respectively. 
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Figure 6 : Lorentz factor ( 7 ) distributions of x?’s that enter (solid lines) or decay in 
(dashed lines) the decay tunnel for (a) Case 1, (b) Case 2 , and (c) Case 3. Hatched 
histograms are for y^’s whose decay photons enter the shower detector. 
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Figure 7: The energy spectra of the deci 
(a) Case 1, (b) Case 2, and (c) Case 3. 
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Figure 8: Scatter plots {ip vs. the arrival time) of detected events for (a) Case 1, 

(b) Case 2, (c) Case 3, and (d) background neutrinos and muons. The resolutions of 
tjj and time measurements are assumed to be up = 0.3° and at = 0.2 ns, respectively. 
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Figure 9: Distributions of the arrival time before (solid lines) and after (dotted lines) the 
bunch convolution for (a) Case 1, (b) Case 2, (c) Case 3, and (d) background neutrinos 
and muons. Arrows indicate the bunch spacing of 25 ns. 
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Figure 10: Explorable range of cr for the Xi decay at the LHC with 300 fb“^. Also 
shown are the predicted curves of cr(x° — 7 G) for ^/F = 3 x 10® to 3 x 10^ GeV (solid 
lines) and cr(x? — 7 a) for Fa/N = lO'^ to 10^° GeV (dashed lines). 







Figure 11: Timing chart of the Xi decay event. The approximate time-of-flight (TOF) of 
Xi, t^oF, is extracted from the time difference between the signal of the shower detector 
and the RF pulse. Note that | OF |/c is the minimum possible TOF on the condition of 

/ 5 < 1 . 
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Figure 12: Distributions of {3^ ^ for (a) Case 1, (b) Case 2, (c) Case 3, and (d) background 
neutrinos and muons. 




Figure 13: Scatter plots of vs. 7 * for (a) Case 1, (b) Case 2 , and (c) Case 3. The 
solid lines correspond to and the dashed lines are for m^o deviated from the 

true value by ± 20 %. 









E/Y. (GeV) 

Figure 14: Distributions of -E^/ 7 * for (a) Case 1, (b) Case 2 , and (c) Case 3. Arrows 
indicate the position of -E^/ 7 * = m^o. The dotted lines show the position of E^/'y^ = 
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Figure 15: Distributions of the off angle tjj for (a) Case 1, (b) Case 2, (c) Case 3, and 
(d) background neutrinos and muons. The cut position for the kinematical analysis is 
also shown. 
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Figure 16: Same as Fig. after applying the cut of -0 > 10°. 
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Figure 17: Distributions of the reconstructed mass for (a) Case 1, (b) Case 2, and 
(c) Case 3. Solid lines show the results of the £t by a gaussian. Arrows indicate the 
position of M q = m^o. 
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Figure 18: Total cross sections for the production of gluinos, squarks and gauginos in pp 
collisions at \/s = 100 TeV as a function of mg. The corresponding values of m^o are also 
shown. Also given in the vertical scale is the corresponding number of events per 1000 
fb-h 


35 


Events per 1000 fb 







cx (m) 



100 200 300 400 500 600 700 800 900 1000 1100 


m-„ (GeV) 

Figure 19 : Explorable range of cr for the Xi decay at VLHC with 3000 fb“^. Also shown 
are the predicted curves of ct{xi 'yG) for ^/F = 3 x 10 ® to 3 x 10 ® GeV (solid lines) 
and cr(x'j’ — ya) for Fa/N = 10 ® to 10 ^^ GeV (dashed lines). 
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Figure 20: The schematic of the shower detector we investigated. It consists of the 
preshower detector (PS), two parts of the electromagnetic calorimeter (EMI and EM2), 
three planes of plastic scintillators (SI, S2 and S3), the lead absorber (LI), and the hadron 
calorimeter (HAD). 
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Figure 21: An example of the simulated electromagnetic shower which is initiated by 
100 GeV photon incident on the shower detector. Solid and dashed lines indicate tracks 
of > 50 MeV electrons and photons, respectively. 
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Figure 22: Distributions of energy depos 
(a) S2, (b) EMI, (c) EM2, and (d) HAD. 
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Figure 23: (a) Distributions of the sum of deposit energy in EMI and EM2 for photons 

incident at 0° on the shower detector with energy E = 20, 50, 100 and 200 GeV. The 
results of the £t by a gaussian are also shown, (b) The sum of deposit energy in EMI 
and EM2 vs. incident photon energy, (c) The energy resolution as a function of incident 
photon energy. The curve is the result of the £t for 0° data. 
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Figure 24: The angular resolution of the reconstructed shower direction as a function of 
incident photon energy. 
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Figure 26: The energy spectra of the prompt muons produced in the decay of bottom 
and charm hadrons for total 2.4 x 10® generated events. The polar angle of the tunnel 
axis 6*axis is taken to be (a) 15°, (b) 20°, (c) 25°, and (d) 30°. 





